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Abstract 

The model proposed by A. Zee (1986) and K. S. Babu (1988) is a simple radiative seesaw model, 



in which tiny neutrino masses are generated at the twodoop level. We investigate a supersym- 
Q-f 

metric extension of the Zee-Babu model under R-parity conservation. The lightest superpartner 

<D : 

—, . particle can then be a dark matter candidate. We find that the neutrino data can be reproduced 

with satisfying current data from lepton flavour violation even in the scenario where not all the 

> 
q^ ■ superpartner particles are heavy. Phenomenology at the Large Hadron Collider is also discussed. 
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Although the standard model (SM) has been successful in describing phenomena below 
100 GeV, the Higgs sector has not been confirmed yet. The Higgs boson is expected to 
be lighter than one TeV from the unitarity argument, so that it can be explored at the 
CERN Large Hadron Collider (LHC). On the other hand, we require new physics beyond 
the SM because of several reasons, such as the quadratic divergence problem, the origin of 
tiny masses of neutrinos, the existence of dark matter (DM), and so on. It is interesting to 
consider a scenario where these problems are simultaneously solved at the TeV scale, as such 
a scenario is directly testable at the LHC or future colliders such as the International Linear 
Collider (LLC) and the Compact Linear Collider (CLIC). In such a case, it is plausible that 
the Higgs sector is closely related to the detail of physics beyond the SM. 

One of the important motivations to consider physics beyond the SM is to explain the 
origin of tiny neutrino masses. The seesaw mechanism is known to be a simple method 



of generating neutrino masses at the tree level, in which tiny masses of the ( 
neutrinos may be obtained from very heavy right-handed neutrinos (type I) 



»E 



eft-handed) 

1 

1[, a heavy 



triplet Higgs boson (type II) |2j], or a heavy triplet fermion (type III)[3j. However, the mass 
scale of these fields is much higher than the TeV scale; naively at around (9(10 6-15 ) GeV, 
unless the coupling constants between lepton doublets and these new heavy fields are taken 
to be unnaturally small. Such a high scale is far from experimental reach. 

Radiative seesaw models, where neutrino masses are generated at the quantum correc- 
tions, are amative attractive scenarios to generate tiny neutrino .nasjfl. Masses of 
new particles in these models can be as low as the TeV scale, so that they are expected to 
be directly testable at current and future collider experiments. One of the characteristic fea- 
tures of these models is an extended Higgs sector. Another feature is the Majorana nature, 
either introducing lepton number violating couplings or introducing right-handed neutrinos. 

The original model for radiative neutrino mass generation was first proposed by A. Zee[4|, 
where neutrino masses are generated at the one- loop level by adding an extra SU(2) L doublet 
scalar field and a charged singlet scalar field with lepton number violating couplings to the 



SM particle entries. Phenomenology of this model has been studied in Ref. lOj. However, 



it turned out that it was difficult to reproduce the current data for neutrino oscillation in 



this original model [111]. Some extensions have been discussed in Ref. [121 ] . 

The simplest successful model today may be the one proposed by A. Zee[5] and 
K. S. Babu[6], in which two kinds of SU(2) L singlet scalar fields are introduced; i.e., a 



singly charged scalar boson and a doubly charged one. These fields carry lepton number 
of two unit. In this model, which we refer to as the Zee-Babu model, the neutrino masses 
are generated at the two-loop level. Phenomenology of this model has been discussed in 



Refs. 13Ml7l. Apart from the Zee-Babu model, there is also another type of radiative see- 



saw models 7-9], where TeV-scale right-handed neutrinos are introduced with the odd charge 
under the discrete Z 2 symmetry. In these models, the Z 2 symmetry protects the decay of 
the lightest Z 2 odd particle, which can be a candidate of DM. This is an advantage of this 



class of models 



I9U2I 



^2l| . On the other hand, in the Zee-Babu model there is no such a discrete 
symmetry and no neutral new particle, so that there is no DM candidate. 

In this Letter, we investigate a supersymmetric extension of the Zee-Babu model. By 
introducing supersymmetry (SUSY), the quadratic divergence in the one-loop correction 
to the mass of the Higgs boson can be eliminated automatically. In addition, a discrete 
symmetry, which is so called the R-parity, is imposed in our model to forbid the term which 
causes the dangerous proton decay. The R-parity also guarantees the stability of the lightest 
super partner particle (LSP) such as the neutralino, which may be identified as a candidate 
of DM. We find that there are allowed parameter regions in which the current neutrino 
oscillation data can be reproduced under the constraint from the lepton flavour violation 
(LFV) data. In addition, this model provides quite interesting phenomenological signals 
in the collider physics; i.e., the existence of singly as well as doubly charged singlet scalar 
bosons and their SUSY partner fermions. Such an allowed parameter region also appears 
even when new particles and their partners are as light as the electroweak scale. We also 
discuss the outline of phenomenology for these particles at the LHC. 

In the original (non-SUSY) Zee-Babu model 5|, la], two kinds of SU(2) L singlet fields uo~ 

(Y = —1) and k (Y = —2) are introduced. The Yukawa interaction and the scalar 

potential are given by 

3 3 

C = -J2 fi/ c u ■ hj^ + - J2 9^Ri e Rj K ~~ ~ V B u~u~K ++ + h.c. -V- V SM , (1) 

where V SM is the Higgs potential of the SM, £ Li are lepton doublets, e Ri are right-handed 
charged leptons, the indices i, j are the flavour indices, the dot product of the fields denotes 
the antisymmetric contract of the SU(2) L indices i.e. I c Li ■ l L j = Yl a B=x e a^ ( Li^Lji an d an 
the scalar couplings with respect to u~ and k, other than u~u~ n ++ are in V . Notice that 
lepton number conservation is broken only by the term of [i B . The neutrino mass matrix is 



generated via two-loop diagrams as shown in Fig. [TJ The induced neutrino mass matrix is 
computed as a 

1 \ 2 3 l6 flB f ik (m e ) k g kl (m e ) l f jl 
TT-j 2^ Z~2 J ( m "' y m e)k\ m ^ Vm e )i\m K ) , (2) 
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where (m e )i are charged lepton masses, and the induced mass matrix (m„)^ is defined in 
the effective Lagrangian as 

3 1 _ 
and I(m u ,m 12 \m 21 ,m 22 \M) is the two-loop integral function defined as 



(3) 



I(m n ,m 12 \m 21 ,m 22 \M) 



M 2 



(p 2 + 77%) (p 2 + m 2 2 ) (q 2 + m 2l ) (q 2 + m 22 ) ((p + q) 2 + M 2 ) 



Following Refs. 22], one can evaluate the function J(77i 11 ,77i 12 |77i 21 , m 22 |M) as 



(4) 



I(m u ,m 12 \m 21 ,m 22 \M) 

M 4 {I(m 12 \m 22 \M) - I(m u \m 22 \M) - I(m 12 \m 21 \M) + I(m u \m 21 \M)} 



(m 2 1 -m 2 2 )(m 2 1 -m 2 2 ) 



where 



J(m 1 |m 2 |M) 
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The function f(x,y) is give by 
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D=y/l-2(x + y) + (x-y¥, 
a + =-(l-x + y + D), r + = -{l + x-y + D), 

a_ =-(1 - x + y-D) , r_ = -(1+X-2/-D) , 



(8) 



Our result for the neutrino mass matrix is consistent with that in Ref. 15[ including the factor 
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FIG. 1: The two-loop diagram relevant to the neutrino mass matrix. 



and Li2(x) is the dilogarithm function defined as 

We note that in the limit of m 12 = m 22 = and m u 



dt 



(9) 
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m^ the above function 
I(m u , m 12 \m 21 ,m 22 \M) has the same form as the function given in Refs. [13l. Il4j. 



I(m u ,Q\m u ,0\m K 
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dx l dy 
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x + (r 
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x + ry 



l)y + y 2 
where r = w? K jm^. One can approximately estimate the above function as 

(r > 1) 

1.98r(r + 0.12) _u ^ , (r < 1) 

Details of the Zee-Babu model have been studied in the literature 



(10) 



I(ra w , 0|m w , 0\m K 



2.8r(r + 0.31) 

,-0.23 



(11) 
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16| . It is known that 



the model can reproduce the present neutrino data with satisfying constraints from the LFV. 
We turn to the SUSY extension of the Zee-Babu model. The SU(2) L singlet chiral su- 
perfields Q C R , fi L , K L , and K C R are added to the superfields in the minimal supersymmetric 
standard model (MSSM), whose details are shown in Table. HI Notice that although the non- 
SUSY Zee-Babu model includes only two SU(2) L singlet scalars these four chiral fields are 
required in the SUSY model. If only fl c R and K L are introduced in the model, their fermion 
components are massless and the model is ruled out. By introducing additional fields Vt L and 
K R such massless fermions can be massive, and furthermore the model becomes anomaly 
free. 



TABLE I: Particle properties of relevant chiral superfields. 
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The superpotential is given by b 



W =^ MSSM + f tl L t ■ LM R + 9ii E?E'K L + \ L K L n c R n c R + X R K R Q L Q L 



+ fki^k^L + \x K K L K c R 



(12) 



where Wmssm is the superpotential in the MSSM. The superfields in the superpotential 
are listed in Table. HJ and the coupling matrices fij and gij are an antisymmetric matrix 
fji = —fij and a symmetric one gji = g^, respectively. It is emphasised that we here impose 
the exact R-parity in order to protect the decay of the LSP, so that the LSP is a candidate 
of the DM. The soft SUSY breaking terms are given by 



^soft — ^MSSM + ^SZB + ^C ' 



(13) 



b Hereafter we omit the summation symbol for simplicity. 



where £mssm represents the corresponding terms in the MSSM, 
Abzb = - M + w R w i? - m -^l ~ M 1 __k' l k l - MI + k* r k r 



+ (- m sfijU*J Li ■ hj ~ m s9ij K Le* m e* Rj - m s \ L K L u R u R - m s \ R K R w L u L 

- B u few> L - B K fi K K L K* R + h.c. ) , (14) 



and 

C c = - CMlh - C d u L ^ u - (Cj^f£ Ll ■ l Lj + h.c. , (15) 

where m s denotes a typical SUSY mass scale, and f^ = — /$• and g^ = g^. £ SZB is the 
standard soft-breaking terms with respect to the new charged singlet fields, u L R and k l r . 



C c contains the terms so-called the "C-terms" 25|, where the scalar component and its 
conjugation are mixed c . 

There are two possibilities in building a SUSY model with the charged singlet fields, 
depending on whether or not the C-terms are switched on in a SUSY breaking scenario d . 
If we assume that C c is absent, tiny neutrino masses are generated only by at least two 
loop diagrams as in the Zee-Babu model. On the other hand, with the term oj^cj^cf)^ tiny 
neutrino masses are dominated by one loop diagrams in Fig. [2] just like in the original Zee 
model[4j. In this Letter we focus on the case where the SUSY breaking mechanism does not 
lead to the soft SUSY breaking C-terms, so that all the neutrino masses are generated at 



the two loop level. The case with the C-term will be discussed elsewhere [27]. 

From the superfields Q C R , Q L , K L and K R , there appear singly charged (Y = —1) and 
doubly charged (Y = —2) singlet scalar bosons, uj r L and n L R , as well as their superpartner 
fermions, namely singly and doubly charged singlinos, Co and k, respectively. The superpo- 
tential and the soft SUSY breaking terms lead to the mass matrix for the singly charged 
scalars in the basis of (u R , uj l ) as , 

} fl ' M% + |/4| - m 2 w tan 2 9 W cos 2/3 B^ n \ 

B u fJ. n M2 + |^|+m 2 v tan 2 ^cos2/3/ 



c The singlet scalar C-terms break SUSY hard, while the terms listed in the C c include non-singlet scalars 

and the quadratic divergence does not occur. 
d Many models derived by N = 1 supergravity do not lead to the C-terms and if they are absent at the cut 



off scale, they do not appear through the radiative corrections 26j. Thus the C-terms are usually ignored 



in the MSSM. On the other hand, it is known that C-terms are induced in some models of SUSY breaking 



such as an intersecting brane model with a flux compactification 24 1. 
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FIG. 2: The one-loop diagram relevant to the neutrino mass matrix with the C-term uj + ^* 



u^ d - 



(y e ) • ■ is the charged lepton Yukawa coupling. 
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and the mass matrix for the doubly charged singlet scalars in the basis of (k l , k r ) as 
., , M 2 _ + \i^ K \ + 2m^ tan 2 9 W cos 2/3 B K /i K 



Mi 



B K fi K M 2 + + \/j, 2 K \ — Im^ tan 2 # w cos 2/3 ; 



(17) 



where tan/3 is a ratio of the two vacuum expectation values of the MSSM Higgs bosons 
as tan/3 = (<p u )/((p d )- As easily seen from the above expressions, uj r and uj l (k l and 
K R ) can mix with each other by the soft-breaking "B-term", {B^ix^)u: r (jJ l ((B k (j, k )k l k* r ). 
The mass eigenvalues of singly and doubly charged singlet scalar bosons are obtained after 
diagonalising their mass matrices M 2 and M 2 by the unitary matrices U u and U K as 
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K K K 



(18) 
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K)2 , 



{m w )il \ 

The mass eigenstates are then given by 

"a = (Ul) al u R + (Ul) a2 uJ L , K a = (Ut) al K L + (Ul) a2 K R , (a = 1, 2) . (19) 

The mass eigenstates of the singlinos are 



UJ, 



UJ 



K 



(20) 



UJ 



Rj 



"R, 



whose mass eigenvalues are given by the SUSY invariant parameters as m^ = /i c and 
m k = n K , respectively. 

The neutrino mass matrix is generated via the two-loop diagrams shown in Fig. [31 which 
can be written as 



1 
where the matrix H kl is a symmetric matrix 



m v)i3 = M^ 2 \2 ^ m ^ H kl{ m e)lfjl , (21) 



H ki = l ^B)abc{UJla{U u) Yi b {U K ) lc g kl I{{'m e ) k , (m w ) a |(m e )„ (m w ) 6 |(m«) c ) 

+ 16^l(C/JL(^)i a {^— /((m g8 ) fc ,m |(m g J I) m ffi |(m K ) a ) 

X X 

+ — k gki I ((me R ) k , m^\(m s )„ m |(mj a ) + g k i—^I{{m s ) k , m a \ (m s ),, m | (mj ) 

r v" v- 

X S — ^fe J ((^e- )fc 5 m <il(^ e )^( m ^)al m s)+fe— Lj ((™ e )fc 5 KJJ(™ e ~ )z,«\>K<) 



m 5 rt 771 c 



(22) 



where the indices a, 6, c run from 1 to 2, the mass eigenstates of the charged singlet scalars, 
(m g )i and (m g ) f are slepton masses, the left-right mixing term in the slepton sector is 
parameterized as (m e ) k X k /m s , I(m n , m 12 \m 21 , m 22 \M) is the loop function given in Eq. (J3J), 
and the other parameters are defined in the relevant Lagrangian as 

C = - 2f ii (U u )ltfP L e j uZ - ^.([/J la e 4 P L e^ a - 2/^P^ - 2f i tfP L u i: e Lj 

~ 2 9i^*R^j P L^ ~ ^L(U K ) la uP L u c K a - 2X L (Uj* la uP L kuj* a - g i j{U^ Xa {m^)fe f Ri e f Lj K a 
- (^ B )a6cW a w 6 K* - m 8 g ij {U K ) la S' Ri ? I iiK a + h.c. , (23) 

with 

(/*b)«*c =A£(l/J ^JlbC^X + MUj 2a (Uj 2b (U K T 2c . (24) 

In the above expression, we assume that there is no flavour mixing in the slepton sector. In 
our model, there are two sources of the LFV processes. One is the slepton mixing which also 
appear in the MSSM. The other is the flavour mixing in the coupling with the charged singlet 
particles. In order to concentrate on the latter contribution to the lepton flavour violating 
phenomena, the usual slepton mixing effect is assumed to be zero. The phenomenological 
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FIG. 3: The contributions to the neutrino mass generations. A type of a diagram (a) is the 
corresponding diagram to the non-SUSY Zee-Babu model. Diagrams (b) and (c) are new type of 
diagram in the SUSY model. 

constraints in our discussion strongly depend on this assumption. If the assumption is 
relaxed, the phenomenological allowed parameters of the model can be changed to some 
extent. Still we think our assumption is valuable to consider in order to obtain some definite 
physics consequences which are relevant to the new particles in our model. 



10 



It is non-trivial whether there is an allowed parameter region in our model except for the 
decoupling limit where masses of all the super partner particles are set to be much larger 
than the electroweak scale. Let us search for the parameter region where the neutrino mixing 
is consistent with the present oscillation data and the LFV constraints are satisfied. 

Flavour violation in couplings between SU(2) L singlet fields and leptons should be large 
in order to generate large off-diagonal elements in the neutrino mass matrix. These large 
flavour violation couplings enhance the LFV processes. In particular doubly charged singlet 
scalar exchange tree level diagram contributes to the ef — > e ^et ef process. The predicted 
decay width of ef — > efe^ef in the model is calculated as[13|, [14| 



r (e+ ->■ eje+e L 



r 1 (m e ) 5 i 
ife 81927T 3 



(U K )* la (U K ) h 



9u9 jk 



[m„ 



(25) 



where Cjf. is a statistical factor as 



a 



jk 



i (j = k) 

2 (J^k) 



(26) 



There can be still large contributions to \x — > e'y, even if the constraint from /i + — > e + e + e~ 
can be avoided. The contribution is from one-loop diagrams. The decay width of e^ — > e^ 
is evaluated as 



with 



Or 



r(e^e j7 ) = ^KX 5 (|^| 2 + |^| 2 ) 



(27) 
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(28) 



(29) 



where (m u )i are neutrino masses, and (mj )j are sneutrino masses. The loop functions Fi(x) 
and F 2 (x) are^ 



Fi(x) 
F 2 (x) 



x 2 — 5x — 2 x In x 

+ 



12(x-l) 3 2(x-l) 4 
2a; 2 + 5a; — 1 x 2 In x 



12(a;-l) 3 2(x-iy 



(30) 
(31) 
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The coupling constants fy only have nonzero values in flavour off-diagonal elements, and 
they tend to be large to reproduce the bi-large mixing. Then the bound from the data 
becomes severe. 

Let us discuss how the LFV processes constrain the parameter space. First of all, the tree 
level diagram contributing to the \i — y eee must be suppressed. The present bound on the 



branching fraction is B(fi + — y e + e + e~) < 1.0 x 10 -12 31], which gives very strong constraint 
on the model parameter space. There are two possible cases to suppress the tree level 
contribution to the /x + — y e + e + e~ . The first possibility is considering heavy doubly charged 
bosons Ki and k 2 . If Q\\ ~ #12 ~ 0.1 is taken, the doubly charged bosons should be heavier 
than 15 TeV to avoid too large contribution. The second option is suppressing a product 
of the couplings \g\29\\\- When the doubly charged bosons are 500 GeV, the upper bound 
on the product |^7i25'ii | is obtained as |^i25 l ii| < 10~ 5 . The contributions to r + —y e + e + e~, 
r + — y e + e + /i~, r + — y [i + fi + e~, t + — y fi + n + fi~, t + — y fi + e + e~, and r + — y e + fi + fi~ can be 
computed in the same manner. These flavour changing tau decays into three leptons are also 
enhanced in the model with tree level contributions. If future tau flavour experiments such 
as the high luminosity B factories[29] would discover a signal of such decays, it could support 
the model. In the phenomenological point of view, the scenario with a light doubly charged 
singlet scalar is attractive because the scenario with such a light exotic particle is testable 
at the LHC. Therefore we have searched for a solution with a suppressed (g^iil and we 
have found that the coupling gu can be taken to be so small that the tree level contribution 
to the /i + —y e + e + e~ process is negligible with reproducing the neutrino oscillation data. In 
such a parameter space, the B(fi + — y e + e + e~) is suppressed by the electromagnetic coupling 
constant compared with B(/j, — y erf), say B{ji + —y e + e + e~) ~ a e B{ji — y ej) where the 



current upper limit is given by B(fi —y cy) < 1.2 x 10~ n 30]. The B(fi + —y e + e + e~) is 
below the experimental upper bound, if the constraint of B(/j, — y ej) is satisfied. 

In our analysis below, we work in the limit of B^fiQ —y and B K fi K —y for simplicity. 
If these terms are switched on, the mixings in the charged singlet scalar mass eigenstates 
take part in the neutrino mass generation. However these mixings do not change our main 
results. In this limit, the mixing matrices U w and U K become the unit matrix, and only u 1 
and k x contribute to the neutrino mass matrix and the LFV. Below we simply write the 
relevant fields as u = oj 1 and k = /%, and their masses are written as m w = (t7\,)i and 
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Following the above strategy, we search for an allowed parameter set. An example of the 
allowed parameter sets is 
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1000 GeV , 




m co 


= 600 GeV 
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= 600 GeV , 


m K = 


= 300 GeV , 



m k = 200 GeV , 

("O2 > m u> , (m K ) 2 > m n ■ (32) 

On this benchmark point, the neutrino masses and mixing angles are given as 

sin 2 0i 2 = 0.33 , sin 2 # 23 = 0.5 , sin 2 13 = 0.0 , 

Am 2 21 = 7.6 x IO" 5 eV 2 , |Am 2 x | = 2.5 x IO" 3 eV 2 , (33) 



which are completely consistent with the present neutrino data: the global data analysis 23] 
of the neutrino oscillation experiments provide sin 2 6*12 = 0.318^o;gjJ, sin 2 823 = 0.50lo;gg, 
sin 2 fl 13 = 0.013j;{j;{jijj, Am 2 x = (7.59+q 1 ; 23 ;) x IO" 5 eV 2 , and (Am^l = (2.40+^ 2 ) x IO" 3 eV 2 . 
Based on this benchmark point, our model predicts B(fi — > erf) = 1.1 x 10" 11 and B(t + — > 
/x + /i + yU") = 1.3 x 10~ 8 , both of which are just below the present experimental bounds. Since 
the LFV is naturally enhanced in the model, the MEG experiment [35j], which is expected 
to achieve B(n — > cj) < IO" 13 in a few years, will cover very wide regions of the parameter 
space. Apart from the bench mark scenario, there can be other parameter sets where the 
neutrino data and LFV data are satisfied. However, we here do not discuss details for such a 



possibility. A more general survey of the parameter regions may be performed elsewhere 27]. 
We turn to discuss collider phenomenology in the model assuming the parameters of 
the benchmark scenario given in Eq. fl32|) . In our model, the new SU(2) L charged singlet 
fields are introduced, which can be accessible at collider experiments such as the LHC 
unless they are too heavy. In particular, the existence of the doubly charged singlet scalar 
boson and its SUSY partner fermion (the doubly charged singlino) provides discriminative 
phenomenological signals. They are produced in pair (k, ++ k or k ++ k ) and each doubly 
charged boson (fermion) can be observed as a same-sign dilepton event, which would be a 
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FIG. 4: Production cross sections of (a) k ++ k and (b) R ++ k , via Drell-Yan processes at 
the LHC (jpp) and the Tevatron (pp). The production cross section at the LHC is evaluated for 
^/s = 7 TeV and y/s = 14 TeV. 

clear signature. In this Letter, we focus on such events including doubly charged particles. 
For the benchmark point given in Eq. ( 132|) . almost all the k decays into the same-sign muon 
pair, k ±±l — >■ A i± /i ± . 

At hadron colliders such as the LHC and the Tevatron, the doubly charged singlet scalar 
k and the doubly charged singlino k are produced dominantly in pair through the Drell- 
Yang processes. The production cross sections for k ++ k and k ++ k are shown as in 
Fig. Ufa) and Fig. 11(b), respectively. The first two plots from above correspond to the cross 
sections at the LHC of y/s = 14 TeV and yfs = 7 TeV, and the lowest one does to that at 
the Tevatron of \fs = 2 TeV. We note that magnitudes of the production cross sections for 
the pair of singly-charged singlet scalars u + u~ and that of singly-charged singlinos u + £j~ 
are (1/4) smaller than those for k ++ k, and k ++ k for the common mass for produced 
particles. The direct search of doubly charged Higgs bosons at the Tevatron gives the lower 
)ound on the mass assuming large branching ratio decaying to muon pairs as m K > 150 GeV 
32J. Such a bound on the mass of doubly charged singlinos is partly discussed in Ref. [33] . 
At the LHC with y^i = 7 TeV with the integrated luminosity £ of 1 fb _1 , about 100 of 
k ++ k pairs can be produced when m~ = 200 GeV, while only a couple of the k, ++ k 
pair is expected for m K = 300 GeV. 

In Fig. [U the distribution of the differential cross section for four muon (plus a missing 
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FIG. 5: The invariant mass distribution of the same-sign dilepton event. The benchmark point 
in Eq. (|32[) is used and the neutralino mass is taken as m^o = 100 GeV. The dashed (red) curve 
corresponds to the events from pp — > k ++ k — > x° K f + X° K i ~~ * X°X°/ i+ / i+ / i ~M~ • The dot- 
dashed (blue) curve shows the contributions from pp — > k ++ k, — > fi + fi + fi~ /j,~ . The solid (black) 
curve denotes total events from the both signal processes. The dotted (green) curve shows the 
background events. For kinematical cut, see the text. 

transverse momentum) final states as a function of the invariant mass M{^i + p + ) of the same- 
sign muon pair is shown assuming the bench mark scenario in Eq. fl32|) at the LHC with 
\/i = 7 TeV. In order to suppress background events, we select the muon events with the 
transverse momentum larger than 20 GeV and the pseudo-rapidity less than 2.5. The signal 
events come from both pp — > k ++ k — > p + p + p~p~ and pp —¥ k ++ k — > X 0f H X° K i~ ~* 
X°X°/ i+ / i+ / i ~fi>~- The M (/i + /x + ) distribution can be a key to explore the phenomena with the 
doubly charged particles. The doubly charged scalar mass and the mass difference between 
the doubly charged singlino and the neutralino are simultaneously determined at the LHC. 
A sharp peak is expected in the M{jj, + jj, + ) distribution at M(/z + /i + ) = m K , because the 
same-sign muon pair from the k decay is not associated with missing particles. On the 
other hand, the doubly charged singlino decays as k — > x° K ~~ ► j(V ^ m the case 
that the lightest R-parity odd particle is a neutralino, x°? which is a DM candidate in the 
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model. In this Letter, we just assume that the LSP neutralino is Bino-like. In our analysis, 
we fix the neutralino mass asm. = 100 GeV. The mass difference between k and x° can 
be measured by looking at a kink at M(yU + /i + ) = m k — m~ in the M(/i + /i + ) distribution. 
The main background comes from four muon events from the SM processes where muons 
are produced via the ZZ, 77 and 7Z production, or a pair production of muons with the 
Z or 7 emission. The expected background is also shown in Fig. |5j The events from signal 
dominate those from the background in the area of M(fi + fi + ) < m~ — m~ and at around 
M(/i + /i + ) ~ m K . The background events have been evaluated by using CalcHEPJ36[]. From 
this rough evaluation, one may expect that the event from the signal can be identified even 
at the LHC with y/s = 7 TeV and C = 1 fb" 1 . As for the case with y/s = 14 TeV, the signal 
to background ratio becomes larger and it will be more promising to explore our model. 

There are other models in which the same-sign dilepton events are predicted. The model 
with the complex triplet scalar fields is an example of such a class of models [34J]. They can 
in principle be distinguished by looking at the decay products from doubly charged fields. In 
our scenario, n ±±: can mainly decay into A i± /i ± , while in the triplet models where the decay 
of doubly charged singlet scalars are directly connected with the neutrino mass matrix, there 
is no solution where only the /j^/j^ mode can be dominant decay mode. The difference in 
such decay pattern can be used to discriminate our model from the triplet models. 

In this Letter, we have not discussed details for DM physics in our model. Assuming 
the Bino-like LSP, we expect that our DM candidate can satisfy the constraints from the 
WMAP data for the DM abundance in a similar way to the case in the similar scenario in 
the MSSM. We still note that the existence of doubly and singly charged particles in our 
model may change the cross section of DM pair annihilation at the one loop level to some 
extent, so that they may affect the DM abundance. The detail is, however, beyond the scope 
of this Letter, which will be discussed elsewhere 27]. 



We also give a comment on the possibility for baryogenesis. There can be several possibil- 



ities to realise baryogenesis in our model, such as using the Affleck-Dine mechanism 37J], low 



energy leptogenesisJ38|], and electroweak baryogenesis (EWBG). In the scenario of EWBG, 



the electroweak phase transition must be of strongly first order. In the MSSM the scenario 



of EWBG turns out to be rather challenging [39|]. On the contrary in our model, such a 
scenario may be natural and realistic. In our model, there are doubly and singly charged 
singlet scalar fields. When they have non-decoupling property J9|, |2l|, |40|], the parameter 
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region of the strong first order phase transition can be much wider than that in the MSSM. 
In addition, there can be many CP violating phases in the model, which are also required 
for successful baryogenesis. 

We have discussed the SUSY extension of the Zee-Babu model under R-parity conserva- 
tion. In the model, it is not necessary to introduce very high energy scale as compared to 
the TeV scale, and the model lies in the reach of the collider experiments and the flavour 
measurements. We have found that the neutrino data can be reproduced with satisfying the 
current bounds from the LFV even in the scenario where not all the superpartner particles 
are heavy. The LSP can be a DM candidate. Phenomenology of doubly charged singlet 
fields has also been discussed at the LHC. 
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